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TECHNICAL  ROTE  2388 

AXISYMMETRIC  SUPERSONIC  FLOW  IN  ROTATING  IMPELLERS 
By  Arthur  W.  Goldstein 


SUMMARY 

General  equations  are  developed  for  isentropic,  frictionless, 
axi symmetric  flow  in  rotating  impellers  with  blade  thickness  taken 
into  account  and  with  blade  forces  eliminated  in  favor  of  the  blade- 
surface  function.  It  is  shown  that  the  total  energy  of  the  gas 
relative  to  the  rotating  coordinate  system  is  dependent  on  the 
stream  function  only,  and  that,  if  the  flow  upstream  of  the  impeller 
is  vortex-free,  a  velocity-potential  function  of  only  the  radial 
and  axial  distances  exists  in  the  impeller. 

The  characteristic  equations  for  supersonic  flow  are  developed 
and  used  to  investigate  flows  in  several  configurations  in  order  to 
ascertain  the  effect  of  variations  of  the  boundary  conditions  on  the 
internal  flow  and  the  work  input.  Conditions  varied  are  prerotation 
of  the  gas,  blade  turning  rate,  gas  velocity  at  the  blade  tips, 
blade  thickness,  and  sweep  of  the  leading  edge. 


INTRODUCTION 

No  method  of  practical  importance  has  been  described  in  the 
literature  for  computing  the  flow  in  the  three-dimensional  region 
between  the  blades  of  the  supersonic  compressor.  This  deficiency  in 
design  technique  has  been  partly  met  by  using  two-dimensional  theory 
to  compute  the  variations  of  the  relative  flow  from  blade  to  blade, 
by  assuming  no  radial  flow  components,  and  by  neglecting  spanwise 
variations.  This  procedure  is  probably  satisfaqtory  for  designs  with 
ratio  of  blade-root  radius  to  blade-tip  radius  close  to  unity,  and 
cylindrical  hub  and  casing  shape.  When  this  ratio  is  small,  the 
spanwise  variations  are  significant.  A  method  is  described  herein, 
which  was  devised  at  the  NACA  Lewis  laboratory,  for  computing  the 
magnitude  of  these  variations,  when  blade-to-blade  variations  are 
neglected,  by  assuming  a  large  number  of  blades  in  the  impeller  and 
dealing  with  the  circumferential  average  flow.  This  technique,  when 
combined  with  estimates  of  blade-to-blade  variations,  is  useful  for, 
designs  with  closely  spaced  blades. 
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The  usefulness  of  this  approach  may  possibly  be  extended  to  super¬ 
sonic  flow  in  a  manner  analogous  to  the  method  described  in  reference  1 
for  subsonic  flow;  the  axisymnetric  solution  Is  first  found  and  then 
used  as  a  base  for  computing  the  blade-to-blade  variations.  The  par¬ 
ticular  method  used  in  the  blade-to-blade  calculation  of  reference  1 
is  series  expansion,  which  is  not  applicable  to  supersonic  flows 
because  of  the  existence  of  solutions  with  discontinuous  derivatives. 

In  reference  2,  the  equations  for  internal  flow  with  zero  blade 
thickness  are  translated  into  the  characteristic  equations  for  the 
derivatives  of  the  stream  function,  and  it  is  shown  there  that  the 
equations  are  hyperbolic  if  the  relative  velocity  is  supersonic, 
provided  the  ratio  of  the  tangential  to  the  axial  velocity  is  a  given 
field  function.  (This  condition  includes  the  special  class  of  blades 
containing  radii  from  the  axis  of  rotation) .  If  the  tangential  velocity 
is  prescribed,  then  the  meridional  velocity  component  is  the  decisive 
factor  in  determining  whether  the  equations  are  hyperbolic. 

The  present  report  treats  the  velocity  components  directly  after 
the  manner  of  reference  3  for  arbitrary  preassigned  blade  shapes.  The 
design  problem  is  solved  by  finding  the  hub  shape  that  will  give  a 
prescribed  velocity  distribution  along  the  casing  (or  any  other 
selected  stream  surface).  In  the  course  of  the  solution,  the  tangen¬ 
tially  averaged  flow  is  completely  determined  in  the  interior  of  the 
impeller . 

The  purpo’se  of  this  paper,  in  addition  to  describing  the  method 
in  a  mathematical  outline,  is  to  examine  the  performance  of  the  tech¬ 
nique  in  computed  examples  and  to  find  if  any  peculiarities  in  the  flow 
arise  from  the  boundary  Conditions  imposed,  and  to  develop  a  background 
of  experience  for  the  design  of  impellers  of  the  type  considered  by  exam¬ 
ination  of  the  effect  of  systematic  variations  of  impeller  design  para¬ 
meters  on  the  internal  flow.  The  parameters  varied  are  -inlet  pre¬ 
rotation,  blade  turning  rate,  gas  velocity  distribution  at  the  casing, 
blade  thickness,  and  sweep  of  the  blade  leading  edge. 

Approximately  160  hours  of  computing  is  required  to  work  out  about 
100  points  for  each  of  the  impellers.  This  time  can  be  reduced  by 
application  of  the  punched-card  computing  technique  to  part  of  the 
procedure. 

DEVELOPMENT  OF  EQUATIONS  FOR  AXIALLY  SYMMETRIC  FLOW 

Equations  describing  the  fluid  flow  are  developed  in  convenient 
form  by  use  of  the  equations  of  motion  and  continuity,  the  conditions 
of  axial  symmetry,  and  isentropic  state  changes.  The  uniformly  distri¬ 
buted  blade  forces,  which  are  inconvenient  in  calculations,  are 
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eliminated  from  the  equations  by  substituting  the  blade  surface 
function,  so  that  the  condition  of  integrability  is  automatically 
satisfied. 


Equation  of  Continuity 


The  continuity  equation  for  steady  flow  is 


V  *(pW)  =  0 

(Symbols  are  defined  in  the  appendix.)  If  axial  symmetry  is  assumed, 
the  equivalent  in  cylindrical  coordinates  is 


3F  (p^)  +  4  (pwr)  = 


o 


The  presence  of  blades  serves  to  modify  these  equations  if  the  space 
filled  by  the  blades  is  taken  into  account.  If  the  impeller  contains 
N  blades  with  thickness  t  measured  in  the  direction  of  rotation, 
the  ratio  of  volume  available  for  fluid  flow  to  the  actual  volume  of 
any  space  element  approaches  (in  the  limit  as  the  number  of  blades 
increases  indefinitely  with  Nt  fixed) 


A  _  Wt_ 

r  "  2jtr 


The  continuity  equation  is  therefore  obtained  by  reducing  the  specific- 
mass-flow  terms  pv  and  pw  by  the  ratio  A/r  to  account  for  the 
area  reduction: 


5 

3? 


(pvA) 
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(pvA)  s=  0 


*  (i) 


Because  none  of  the  scalar  quantities  is  a  function  of  6>  equa¬ 
tion  (l)  implies  the  existence  of  a  stream  function  ^(r^z)  such  that 

=  PWA.  1 


=  -  PVA 


> 


(2) 


*e  =  0 


(subscripts  indicate  partial  derivatives)  or 
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7^  =  pAj  X  V  =  pAj  X  W 


(2a) 


Equation  of  Motion 

Substitution  of  energy  for  pressure.  -  The  equation  of  motion  for 
steady,  nonviscous,  isentropic  flow  in  terms  of  the  relative  velocity  is 

2 

V  T-WX  VXW+-7p  =  0)2R+2W  xH+F 

^  p 

p 

where  o>  R  is  the  centrifugal  force  and  2W  X  Cl  the  Coriolis  force. 

For  isentropic  conditions,  the  substitution 


iyp  =yh 


(3) 


(4) 

is  a  measure  of  the  total  energy  of  the  gas  particles  in  the  relative 
system  of  coordinates  and  remains  constant  (as  will  be  shown  for 
frictionless  flow)  for  any  particle;  energy  is  interchanged  among  the 
forms  of  enthalpy,  kinetic  energy,  and  potential  energy  (  -  o?r^/2) 
resulting  from  the  centrifugal-force  field.  In  reference  2,  E  is 
shown  to  be  constant  for  any  particle  for  frictionless  flow  and  no 
heat  transfer,  with  a  finite  number  of  blades  in  the  impeller. 


may  be  made,  so  that  the  equation  is  then 
Tr2  2  2 

,  ,  W  (D  r  ,  . 

Vlh  +  —  -  ~2~  y=wx  (7XW  +  2Q)  +  F  =  wxy  xv  +  f 
The  quantity 


„2  2  2  2 

E  =  h  +  -p-  -  —  tj"'  =  h  +  -x-  -  cor(u  +  oar) 


In  order  to  show  that  E  is  constant  for  any  particle  in  axially 
symmetric  flow,  the  direction  of  the  fielcf  force  is  first  identified. 

If  the  blades  are  considered  to  increase  indefinitely  in  number  and  to 
retain  their  camber-line  shape,  the  forces  exerted  by  the  blades  are 
merely  the  pressure  forces  exerted  normal  to  the  blade  surfaces,  which 
in  the  limit  approach  uniformly  distributed  forces  in  the  body  of  the 
fluid.  The  relative  flow  lines  are  always  contained  in  one  of  the 
family  of  blade  surfaces,  so  that  the  field  or  blade  forces  must  be 
normal  to  the  velocity.  If  the  scalar  product  of  equation  (3)  is  taken 
with  the  relative  velocity,  the  result  is  then 
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W  *  VE  =  0 

which  indicates  that  E  does  not  change  along  any  flow  line. 

Elimination  of  hlade  force .  -  The  "blade  force  can  he  eliminated 
from  the  equation  of  motion  if  the  equation  for  the  hlade  surfaces  is 
assumed  known.  The  force  may  then  he  evaluated  after  the  velocity  com¬ 
ponents  are  known. 

Under  the  condition  of  axial  symmetry,  the  equation  for  the  family 
of  hlade  surfaces 

a(r,  0,  z)  =  constant 

may  he  simplified.  If  the  constant  is  changed,  there  is  a  shift  from 
one  hlade  to  another.  Because  of  axial  symmetry,  this  shift  is  equi¬ 
valent  to  rotating  the  original  hlade  by  a  fixed  angle.  Therefore,  the 
equation  for  the  surfaces  can  he  put  into  the  form 

a  =  6  +  f (r,z)  (5) 

Because  the  relative- velocity  flow  lines  are  .contained  in  the  hlade 
surfaces,  W  is  normal  to  y cl  as  well  as  normal  to  yi|r.  Therefore, 
the  relation  among  the  three  vectors  can  he  shown  to  he 


W  =  X  ya 

From  the  condition  of  axial  symmetry. 


(6) 


c)E 

3e  = 


r  j  *y  E  =  0 


Because  both  yE  and  y\|r  are  normal  to  j  and  W,  they  must  he 
parallel,  and  consequently  E  is  known  to  he  a  function  of  ilr  only 
(reference  4);  that  is, 


E  =  E(t) 

VE  =  Vt  re¬ 
insertion  of  equations  (4),  (6),  and  (7)  into  equation  (3)  yields 
(•without  resolution  of  the  term  V  x  V) 


(?) 
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dE  ,  /  r 

5j?vt  =  (^v» 


V 


*4  4 


Va  •  V  X  Vjvt  +  F 


The  blade  pressure  forces  F  are  normal  to  the  blades  and  therefore 
parallel  to  va)  this  relation  thus  includes  only  two  nonparallel 
vectors  Va  and  y'lb  so  that  the  components  must  separately  satisfy 
the  relations 


V\(r  •  V 


X  ^va-  = 


(r  j  •  W  X  v  x  V)ya 


(8) 


and 


dE 

dxjf 


4  W  •  V  XV 

pA 


(9) 


Special  conditions.  -  Two  special  conditions  of  technical  import¬ 
ance  are  considered:  In  high-speed  machines,  blades  with  radial  ele¬ 
ments  are  desirable  in  order  to  eliminate  bending  loads  resulting  from 
centrifugal  forces.  Another  frequently  occurring  circumstance  is  that 
of  zero  vorticity  upstream  of  the  impeller. 

For  blades  with  radial  elements,  the  vector  i  lies  in  the  blade 
surface  and  therefore 

i  •  ya  =  f  r  =  0 

Therefore  f  is  a  function  of  z  alone.  From  equation  (6), 


0  =  W  •  Va  =  -  +  vfr  +  wfz 


(10) 


so  that  for  radial  “blade  elements 


u  = 


w  r 


df 

dz 


(10a) 


This  relation  provides  considerable  simplification  of  the  equations 
developed  later  for  solving  the  flow  problem. 

For  vortex- free  flow  upstream  of  the  impeller ,  equation  (3)  is 
reduced  for  conditions  there  by  setting  Q=F=^X  V  =  0  and  substi¬ 
tuting  V  for  W  to  get 
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h  +  —  =  constant 

Also,  zero  vorticity  gives,  for  axial  symmetry, 

0  =  7  X  V  =  -  !  ^  (ii  +  oxr)  +  j  ^  -  3^  +  -  ^  (ur  +  cor  ) 

From  the  first  and  last  terms  and  axial  symmetry, 

p 

ur  +  our  =  constant 

► 

Then  the  value  of  the  total  energy  E  relative  to  the  impeller  is 

^  t_  ,  W2  oo2r2  .  V2  t  2\ 

E  =  h  +  -g —  — g —  =  h  +  -g —  oo(ru  +  oar  ) 

which  is  constant  upstream  of  the  impeller.  Then  hy  equation  (7)  E 
remains  constant  in  the  impeller,  so  that  the  equation  of  motion  (9) 
becomes 

7  a  •  V  X  V  =  0  (i: 

If  the  system  of  coordinates  used  is  r,  z,  and  a  for  expressing 
equation  (ll) ,  there  results  (reference  4) 

A  (rV  -7a  X  k)  =  A  (rV  •  i  X  7a) 

02  r,a  or  z,a 

The  subscripts  indicate  which  variables  are  held  constant.  This  equa¬ 
tion  and  the  fact  that  the  quantities  in  parentheses  are  functions  of 

r  and  z  only  establish  the  existence  of  a  potential  function  of 
r  and  z  only  ($(r,z))  such  that 

=  rV.-VaX  k 

z,a 


=  rV •  i  XV  a 


Sa 


0 
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The  last  relation  indicates  that  with  r  and  z  constant,  ®  is 
constant,  and  this  is  true  regardless  of  the  third  independent 
variable,  which  may  be  taken  as  0.  The  three  equations  for  $  are 
then  put  into  vector  form: 

V  *  =  i(rV  -ycc  X  k)  +  k(rV  *  i  X  ya)  (12) 

The  absolute  velocity  may  then  be  shown  to  be 

V  =  y$>  +  r(u  +  as r)y<x  (13) 

and  from 

V  *y<x  =  W  *ya  +  carj*ya  =  as 
equation  (13)  gives  for  the  moment  of  momentum  p, 

p  =r(u  +  osr)  =  m  ~ V  ^  7a  (14) 

(yar 

so  that 

y,  '  (13a) 

(V») 

The  potential  $  can  be  used  as  an  alternative  to  the  stream  function 
in  solving  for  the  flow.  Equation  (13)  indicates  that  the  projection 
of  V$  on  the  blade  surface  is  equal  to  the  projection  of  the  absolute 
velocity  on  the  blade  surface.  The  equation  for  \|r  and  that  for  $ 
are  obtained  in  the  following  manner:  The  velocity  V  is  eliminated 
from  equation  (9)  by  means  of  equation  (6)  and  there  is  obtained 


fv  \  0  _  pA  dE 

va  #  v  x  \£R  v*  x  va )  +  =  -  ~  d>F 

which  is  also  applicable  in  the  general  case  with  nonzero  vorticity 
upstream  of  the  impeller. 
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In  order  to  obtain  an  equation  for  $,  equation  (13a)  is  substituted 
into  equation  (2a) : 

Vt  =  pArye  x  v  =  pArye  xy$  +  pApiye  xya 

Then 

0  =  yxy\|r  =  yx  jpAry©  x  ( y  $  +  pya)] 

Although  obtaining  a  solution  in  terms  of  \|r  is  feasible  whether 
dE/dty  is  zero  or  not,  when  dE/df  =  0,  $  may  prove  more  convenient. 
In  regions  where  W^/a^  is  nearly  1.0,  because  |W|  is  a  double  - 
valued  function  of  the  stream- function  derivative  p  |w|  ,  it  will  be 
uncertain  whether  the  flow  is  subsonic  or  supersonic.  A  solution 
obtained  from  a  network  would  require  a  very  fine  spacing  in  regions 
where  |w|  Z  a. 


General  Equations  in  Cylindrical  Coordinates 

In  general,  where  f r  ^  0  and  dE/di|r  /  0,  the  equations  for  the 
internal  flow  are  solved  simultaneously  for  the  two  velocity  components 
v  and  w,  when  u  has  been  eliminated  from  the  equations  by  using 
equation  (10).  The  equations  are  translated  from  their  vector  form 
to  scalar  notation  by  utilizing  cylindrical  coordinates.  Then  a 
further  transformation  is  employed  with  characteristic  quantities  as 
the  independent  variables. 

The  continuity  equation  (1)  is  first  written  in  vector  form  (the 
condition  of  axial  symmetry  being  used)  and  multiplied  by  r/pA: 

^V*(prW)  =  v,W  +  W,?  +  IW*  V(A/r)  =  0 


For  iserrfcropic  flov 
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Introduction  of  E  t>y  equation  (4)  and  use  of  W  *yE  =  0  produces 


This  expression  finally  reduces  to 


(*  "  a2)  ”z  "  I2  ^z  +  “r*  +  f  "  a2)  Vr  ’  a2  (2  )r '  a2  (2); 


5(lt!^)  +  i[T^l  +  ¥^1-0 


In  cylindrical  coordinates,  equation  (9)  becomes 


^  ^  .  dE 

wr  -  vz  +  fr  -  f z  ^r  =  ^  d* 


(15) 


(16) 


If  equation  (14)  is  now  substituted  into  equations  (15)  and  (16) ,  there 
results  for  continuity 


w2  ruw  fz 

11  "  12  +  I  Vz  + 


ruv  fz  w  ]  I  ruw  fr  w  ^ 

a2  '  ¥r  +  V  a2  '  a ^  Vz 


v2  ru  frv 
1'?  +  — 


vr  +  ^  (w2  fzz  +  2wv  frz  +  V2  frr  I  + 


-) 


2  2 


-  +  I  [v  S  »/r)  +  w  4  <6/r)]  -  0 


and  for  the  equation  of  motion 

^1  +  r  fz  ^  wr  +  r  fzfr  (vp  -  wz)  -  (l  +  r  fr  )  vz  + 

f 


r2f  2 
z 


a  /fr\  a  /fr 
dr[fJ  ^[fz 


-  2  tx  -  pA  0 

r  ^  d\|r 


(17) 


(18) 


« 
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Use  of  Characteristic  Variables 

Statement  of  basic  equations.  -  The  method  of  reference  3  is 
applied  to  equations  (17)  and  (18) .  The  object  is  to  find  two  alter¬ 
native  equations  with  new  variables  \  and  tj,  instead  of  z  and  r, 
such  that  each  equation  (corresponding  to  equations  (17)  and  (18)) 
contains  derivatives  with  respect  to  only  one  of  the  independent 
variables .  Two  additional  equations  are  required  to  find  the  position 
coordinates  z  and  r  in  terms  of  the  characteristic  coordinates 
|  and  r] .  There  are  two  characteristic  slopes  dr/dz,  £+  and 

the  first  corresponding  to  \  varying  and  T)  constant,  the  second 
to  £  constant  and  T|  varying.  If  £  represents  either  of  these 
slopes  and  cr  represents  either  of  the  coordinates  corresponding  to  the 
appropriate  the  results  are 


r0  =  t  z0 

-  r2  frfz  -  ~  r2  (va)2±Vr2  (qa)2  (m2  -  1) 

^  EL 

(the  upper  sign  of  ±  refers  to  the  lower  to  £_)  and 


^  B 


Jw0  +  Gv0  +  Lza  =  0 


where 


.2,  2  v2  _2 


B  =  1  +  rcf/  -  ^  r2  foa)2 


M 2  _  (u  +J__  +  vf) 
a2 

(w  -  rufz)2 

J  — - 5 - (1  +  r2f  2) 


s  (M2  -  1)  r2frfz  +  22  r2  (ya)2±  ^r2  (ya)2  (M2  -  1) 


1  =-  ( 
B,  -3  f 


.2*  2 


r*  f  f  t  +  1  +  r^f  “ID.  + 
1  r  z  z 


h 


-2  (rufz-w) 
a 


•T') 


2f  __  +  +  v2f 

zz  rz  rr 


\  ^  u2  I  .  vr  d(A/r)  ,  wr  d(A/r) 

)  x-^-+x^ 


(19) 

(20) 

(21) 


D2  =  r2f z2 


dE 

2  (u  +  mr)fz  -  pA  ^ 


>(21a) 
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Equation  (20)  for  £  shows  that  the  system  is  hyperbolic  if  M2  >1. 

In  the  homogeneous  equations  (D^  =  Dg  -  0) ,  the  coefficients 
L+  and  L_  (corresponding  to  £+  and  £_)  are  zeroj  consequently, 
in  the  hodograph  plane  the  slope  of  the  characteristics  would  he 

T 

given  hy  ^  =  -  ■g.  In  plane  flow,  j/G  is  a  function  of  v  and  w 

only,  and  this  equation  may  be  integrated  and  the  hodograph  solu¬ 
tions  tabulated  in  advance  of  application  to  any  particular  problem 
that  may  arise.  In  the  present  circumstances,  however,  j/G  also 
involves  the  position  coordinates ,  and  D-^  4  o,  d2  4  0.  A  corre¬ 
spondingly  easy  method  of  solution  hy  simple  waves  is  therefore 
unavailable.  Instead,  the  four  (a  =  |  or  T])  equations  (19)  and  (21) 
are  integrated  simultaneously  in  a  step-by-step  procedure. 

Method  of  computation.  -  In  figure  1,  conditions  are ‘ supposedly 
known  at  P^  and  P^,  and  conditions  are  to  be  computed  at  P^.  Then 

equation  (19)  in  difference  form  is 

r3  ”  ri  =  (z3  -  zj)i+  (22) 

r3  ”  r2  “  <z3  •  Z2^-  (22a) 

This  system  may  be  solved  for  r3  and  z^.  Then  similarly  for  equa-' 
tion  (21),  the  difference  form  is 

j(w3  -  w^)  +  G+  (vj  -  v^)  +  L+  (z3  -  z^)  =  0  (23) 

j(w3  -  w2)  +  G_  (v3  -  vg)  +  L_  (z3  -  zg)  =  0  (23a) 

The  subscript  +  indicates  that  q  is  constant  and  the  subscript  - 
indicates  that  £  is  constant.  In  equation  (22)  and  (23a)  the 
coefficients  £,  J,  G,  L  are  averages  between  values  for  the  points 
P-j-  and  P3.  For  the  first  approximation,  only  the  values  at  P^  are 

known  and  these  values  are  used. 

When  approximate  velocity  components  and  position  coordinates  are 
thus  determined  for  .  P3,  the  coefficients  £,  J,  G,  L  may  be  computed 

there  and  averaged  for  a  better  approximation  of  the  position  coordi¬ 
nates  and  the  velocity  components  at  P^.  This  process  is  repeated  as 


2137 


MCA  TN  2388 


13 


often  as  necessary  to  arrive  at  an  answer  of  satisfactory  accuracy.  The 
smaller  the  interval,  the  quicker  the  convergence  to  the  solution.  In 
each  stage  of  the  process,  if  dE/di|r  ^  0,  then  \|f,  must  be  computed 
from  ,J 


r~ 

^3  -  t  =  I  +  ^z)dz  =  (tyr  + 

From  the  value  of  i|r5,  dE/d\|r  at  P3  may  be  found. 


BOUNDARY  CONDITIONS 
General 


The  discussion  in  the  previous  section  implies  that  boundary  con¬ 
ditions  are  known  on  a  curve  C  containing  points  such  as  P  and  P  , 

J-  Ci 

and  that  the  solution  is  then  extended  into  a  region  containing  points 
such  as  P^.  The  shape  of  the  blades  and  E(i|r)  are  assumed  known  in 

that  region.  On  the  curve  C,  there  are  known  in  addition  the  values 
for  gw  and  v,  from  which  all  other  variables  (such  as  dE/d\jr,  i|r,  u, 
p,  a  ,  and  so  forth)  may  be  computed.  Suppose  the  initial  curve  is 
such  that  the  characteristic  curves  intersect  it  at  two  points  such  as 
A  and  D  or  A  and  E  in  figure  2.  Then  by  equation  (21),  there  is  a 
relation  between  the  variables  at  A  and  D  that  does  not  permit  the  inde¬ 
pendent  specification  of  the  variables  v  and  w  at  A  and  D.  Simi¬ 
larly  if  two  characteristics  passing  through  two  boundary  points  reach 
a  third  boundary  point,  then  another  degree  of  freedom  is  lost.  In 
figure  2,  if  two  variables  are  prescribed  in  AB  and  EF,  then  one  variable 
may  be  prescribed  on  BC  and  DE,  and  none  on  CD. 


Tip  Conditions 

The  supersonic  impeller  has  a  high  pressure  ratio  per  stage  and  a 
large  air-handling  capacity  per  unit  of  frontal  area.  In  order  to  main¬ 
tain  the  large  capacity,  the  tip  radius  at  the  entrance  should  be  at 
least  as  large  as  the  exit.  In  order  to  provide  a  maximum  work  output 
with  a  minimum  turning,  the  exit  tip  radius  should  be  at  least  as  large 
as  the  entrance  tip  radius.  A  reasonable  shape  for  the  tip  streamline 
is  therefore  a  surface  of  constant  radius,  although  smaller  radii  between 
the  entrance  and  the  exit  might  be  useful  in  order  to  affect  the  internal 
flow.  If  the  casing  shape  is  selected,  then,  as  it  is  also  a  streamline, 
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the  ratio  v/w  is  specified.  Selection  of  the  velocity  or  a  velocity 
component  remains.  This  selection  may  be  made  by  consideration  of 
the  blade- to-blade  velocity  variations  and  by  using  some  approximate 
method  such  as  that  outlined  in  the  following  paragraphs. 

In  three-dimensional  motion  (finite  number  of  blades),  Kelvin's 
theorem  states  that  if  the  absolute  velocity  is  irrotational  upstream 
of  the  impeller,  it  will  remain  so  inside  the  impeller;  that  is, 

V  X  V'  =  V  X  W  +  2oik  =  0 

where  the  prime  indicates  velocities  of  the  nonsymmetric  flow.  Then 


V\|r  •  VX  V'  =  0 


may  be  reduced  to 


3 


(W  •  q)  =  r  +  (u'+2oir)  q  •  i 


(24) 


where  chi'/ds  _is  the  arc-length  partial  derivative  of  u' 
direction  of  q.  This  equation  is  now  integrated  from  the 
surface  of  one  “blade  to  the  suction  surface  of  the  next  in 
r  =  constant >  z  =  constant.  The  variable  6  increases  by 
2rtA/Nr,  so  that 


along  the 
pressure 
the  direction 
the  amount 


(V  •  q)  -  (w  .  q) 

s  p 


2jtA 

u’d0  +  2mrq  •  i 


If  it  is  now  assumed  that  the  values  of  ch^/ds  and  of  u'  averaged 

across  the  channel  are  equal  to  the  values  obtained  from  the  axisymmetric 

solution,  and  further  that  (W* •  q)  =  |w T  |,  then 

s  *  s 

H  -  l%l  -  [& ♦  (? H  if,]  ¥  (“*) 

if  the  relation 

v/|w|  =  i  •  q 
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is  utilized.  In  this  equation  all  the  velocity  components  without  sub¬ 
scripts  are  obtained  from  the . axisymmetric  solution.  If 

|W'S|  +  |Wp|  =  2  |w| 

is  assumed  or  some  other  connection  between  |W'g|,  |w  | ,  |w|  ,  and 

the  impeller  geometry,  then  the  blade  surface  velocities  at  the  casing 
may  be  obtained  with  that  relation  and  equation  (24a) .  The  distribu-  ’ 
tion  of  |W|  at  the  casing  may  therefore  be  selected  on  the  basis  of 
obtaining  a  desired  blade-surface  velocity  distribution  or  distribution 
of  u. 


Suppose  that,  in  some  manner,  all  the  variables  have  been  selected 
at  the  case  (curve  BC  in  fig.  3(a)).  The  flow  is  then  uniquely  deter¬ 
mined  in  the  region  BP^C  enclosed  by  the  casing  and  two  characteristic 

lines.  If  the  flow  is  also  prescribed  at  the  leading-edge  surface  BB', 
then  the  flow  is  determined  in  the  extended  region  B'BCP, P-.  When  the 

La 

characteristic  line  BP3  falls  inside  this  region  (fig.  3(a)),  no  diffi¬ 
culty  results.  If,  however,  the  characteristic  B'P4  intersects  the 

casing  inside  the  blade  region  (fig.  3(b)),  then  the  velocity  may  not 
be  prescribed  on  the  curve  BP4  (refer  to  discussion  of  fig.  2),  or 

alternatively,  only  one  of  the  velocity  components  may  be  prescribed 
on  BB'.  (BB'  is  not  a  streamline.)  Similar  considerations  are  valid  on 
the  curve  CC'  at  the  exit. 

The  selection  of  values  for  the  variables  v  -w  on  the 

entrance  edge  BB'  must  be  made  by  joining  the  solution  for  the  flow  in 
the  region  ABB 'A'  with  that  for  the  flow  in  the  blade  region  BCC'B'  . 
Because  the  blades  may  deflect  the  air  at  BB'  discontinuously  and 
because  the  effect  of  thickness  may  be  felt  by  the  air  in  a  discontinu¬ 
ous  manner,  the  discontinuities  in  the  flow  should  be  evaluated.  A 
method  for  making  this  estimate  follows: 

The  angle  p  ‘is  defined  as  the  angle  between  the  meridional 
projection  of  the  blade  leading  edge  and  the  z-axis  (fig.  3(a)).  The 
conditions  outside  the  blade  region  are  designated  by  subscript  zero, 
inside  the  blade  region  by  no  subscript.  Then  across  the  surface  BB' 
the  equation  of  continuity  yields 

po  <wo  sin  P  "  vo  cos  P)  =  f1  "  P  (w  sin  P  -  v  cos  p) 


(25) 
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The  discontinuity  in  blade  thickness  can  cause  no  force  in  the  entrance 
surface  BB',  and  the  blades  can  exert  no  force  in  the  blade  surface. 
Therefore,  no  force  and  no  discontinuity  in  momentum  exist  in  the 
direction  of  the  curves  representing  the  intersection  of  these  two 
surfaces;  that  is 


(W  -  W0)  *(i  cos  p  -  k  sin  p)  x  700  =  0 


or 


r(u  -  uQ)  (fr  sin  p  +  cos  p)  =  (v  -  vQ)  sin  p  +  (w  -  wQ)  cos  p 
There  is  also  available 

u  =  -  r  (vfz  +  vfr) 

If  the  process  is  assumed  isentropic,  then 


i  -  te) 


/  2  21 

^.(1  +  !2_i 1 

\X  2h„ 


r-i 


(26) 

(27) 


(28) 


is  obtained  by  utilizing  the  fact  that 


„  .  W02  ^  w2  ca2r2 

E-hQ+2"*  2  2 


and  that 


There  are  then  available  equations  (25),  (26),  (27),  and  (28)  for 
finding  p,  u,  v,  and  w.  The  initial  values  of  fr  and  fz  may  he 

chosen  from  these  equations  to  obtain  desired  relations  between  the 
velocity  components  outside  and  inside  the  blade  regions. 


Applications  of  Method  of  Characteristics 

In  the  examples  computed  in  this  section,  the  conditions  for  which 
solutions  are  obtained  vary.  The  purpose  of  these  solutions  is  to  find 
the  effect  of  some  variations  of  prescribed  conditions  on  the  resultant 
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flow.  For  all  the  examples;  dE/d\|r  =  0;  the  flow  is  vortex- free 
upstream  of  the  impeller;  at  the  leading  and  trailing  edges  the  radial 
component  v  of  the  velocity  is  zero;  and  the  blades  contain  radial 
elements  (fr  =  0)  .  Relative  Mach  number  at  the  conjunction  of  the 

leading  edge  of  the  blade  and  the  Impeller  casing  is  1.94.  Except  in 
example  VI;  the  value  of  the  blade  angle  with  respect  to  the  axial 
direction  is  zero  at  the  exit  (f  =  0)  .  The  other  conditions  are 
summarized  in  the  following  table: 


Entrance 

Example 

Pre¬ 

rotation; 

(u+ar)r 

Axial 

inlet 

velocity 

Tip 

speed 

Blade 

angle 

at 

tip 

(deg) 

Case 

velocity 

Turning 

rate 

Sweep- 

back 

(deg) 

Blade 

thickness 

I 

0 

0.85a-^ 

1.68a^  | 

-61.8 

(a) 

0  to 
maximum 
to  0 

0 

0 

II 

-0.3a1rt 

0. 835a^ 

1.44a1 

-64.4 

1.94a1 

0  to 
maximum 
to  0 

45 

0 

III 

-0.3a  r 

_L  Z 

0.835a^ 

1.44ai 

-64.4 

1.94ai 

Twice 

example 

II 

45 

0 

IV 

-°-3airt 

0. 83 5a^ 

1.44a' 

0>) 

Twice 

example 

II 

45 

0 

V 

-0.3a  r 

X.  z 

;  0.835a^ 

1.44a_L 

(b) 

Twice 

example 

II 

45 

(c) 

VI 

0. 83 5a^ 

1.44a  ! 

-64.4 

(<*) 

0  to 
maximum 

^10 

0 

s. 

Determined  by  conditions  at  leading  edge  to  z/r^  =  0.356;  held  con¬ 
stant  at  2.1a  from  z/r^  -  0.356  to  1.0. 

^Constant  at  1.94a  to  z/r^  =  0.25;  decreased  to  1.4a  1.6a^)  at  exit. 

cZero  from  z  =  0  to  z/r^  =  0.25;  from  z/r^  =  0.25  to  z/r^  =  0.5 
the  correction  for  blade  thickness  is 


A  =  r  exp 


0.1 

r 


1  -  cos  8 jt 


(t  -  °-25)l' 


^Determined  by  leading  edge;  initial  value;  1.94a. 
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Example  I.  -  The  results  of  the  computation  for  the  first  example 
are  suj5narized”in  figure  4,  which  shows  streamlines,  contours  of  con¬ 
stant  Mach  number  and  constant  energy  level,  and  some  of  the  character¬ 
istics  of  the  £ -family,  which  intersect  the  leading  edge.  Because  of 
the  intersection  of  the  characteristics  from  the  leading  edge  with  the 
tip  of  the  blades,  there  is  no  control  over  the  Mach  number,  as  the 
blade  shape  was  prescribed  by 


d9 

dz 


=  -  r 


-  ! 
dz 


04 


(i  ♦  cos  *  i) 


The  Mach  number  at  the  case  rises  to  a  maximum  of  2.1  at  z/r^  =  0.356, 
and  is  then  assigned  a  constant  value.  By  equation  (4),  the  stagnation 
enthalpy  H  (=h  +  V2/2),  which  measures  the  energy  of  the  gas, 
increases  with  z  at  the  rate 


SH 

= 


cor 


du 


At  the  case 


v  =  0,  and  therefore,  from  equation  (14)  and  W' 


u 


,  2 
+  w  , 


dH  __  co r 
^  ”  |V| 


rw3  d2f\ 
l«l  dz2/ 


This  expression  shows  that  the  rate  of  work  input  is  increased  by 
increasing  the  rate  of  blade  turning  and  decreased  by  increase  in  W. 
Figure  4  indicates  that  the  rate  of  turning  of  the  blade  is  too  small 
in  the  beginning  of  the  impeller  and  that  as  a  consequence  nearly  all 
the  work  input  is  concentrated  in  the  back  end  of  the  impeller.  This 
condition  is  aggravated  by  the  acceleration  of  the  gas  in  the  initial 
portion  (u|Wlz<0  because  u<  0)  .  Increasing  the  turning  rate, 
however,  will  open  the  flow  area  and  cause  |W|  to  increase  unless  the 
area  is  controlled  by  blade  thickness. 

A  great  disparity  of  work  exists  between  the  root  and  the  tip;  the 
values  for  ATT j are  0.65  and  2.5  at  the  exit.  For  isentropic  com¬ 
pression,  the  corresponding  ratios  of  outlet  to  inlet  stagnation 
pressures  are  2.48  and  9.0  for  the  root  and  the  tip,  respectively.  The 
Mach  number  at  the  root  at  the  exit  (z  =  r^,  r  =  0.507  r^)  is  in  the 

neighborhood  of  2.25,  and  a  maximum  value  of  about  2.5  is  reached  at 
z  =s  r^_,  r  a  0.70  r^.  Subsonic  Mach  numbers  are  reached  in  the  region 

0  <  z/r_|_  <0.3.  The  solution  could  not  be  extended  into  this  region. 
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If  the  velocity  were  caused  to  decrease  at  the  exit  (z  =  r^)  by 
decreasing  the  velocity  at  the  case  (r  =  r^)  in  the  rear  part  of  the 
impeller  (0.356 < z/r^ < 1.0) ,  then  the  mass-flow  density  at  the  exit 
would  increase  with  the  decrease  in  velocity  and  the  required  flow  area 
would  decrease .  Termination  of  the  hub  would  be  expected,  not  at 
r  =  0.507  r^.,  but  at  some  larger  value  where  the  work  output  would  be 
higher.  Another  method  of  reducing  the  blade  height  at  the  exit  and 
thereby  reducing  the  work  variation  there  would  consist  in  relaxing 
the  condition  of  v  =  0  and  utilizing  some  other  relations  that  would 
correspond  to  lower  Mach  numbers.  Two  variables  or  relations  could  also 
be  prescribed  at  the  exit,  but  conditions  on  the  blade  tip  (r  =  r^) 

could  not  be  specified  for  some  region  near  the  back  of  the  impeller. 

In  this  respect  the  entrance  and  exit  conditions  would  be  similar. 

Example  II:  Prerotation  of  gas.  -  In  order  to  maintain  supersonic 
flows  throughout  the  impellers  so  that  the  flow  could  be  computed,  the 
relative  entrance  Mach  number  at  the  blade  root  for  example  II  was 
increased  by  assuming  an  initial  counterrotation.  The  rotative  speed 
had  to  be  reduced  from  1.68a^  to  1.44a^  to  maintain  the  work  output, 

the  relative  Mach  number  at  the  leading-edge  tip,  approximately  the  same 
blade  angle  there,  and  the  same  axial  velocity  component  at  the  inlet. 
This  reduction  had  a  further  beneficial  result  of  reducing  the  absolute 
Mach  number  at  the  exit  for  a  given  relative  velocity  axially  directed. 
The  condition  of  z  =  constant  for  the  leading  edge  cannot  be  main¬ 
tained  for  shockless  entry  with  radial  blade  elements  and  a  counter¬ 
rotating  free  vortex  entering  the  blade  system.  The  sweep  angle  p 
chosen  for  this  example  was  45°,  enough  to  permit  free  selection  of 
W  =  1.94a  over  the  entire  blade  tip,  because  the  characteristic  curves 
from  the  leading  edge  do  not  interfere  with  boundary  conditions  on  the 
casing.  The  reduction  in  velocity  at  the  casing  reduces  the  blade 
height  at  the  exit  and  the  work-output  variation  there. 

In  order  to  examine  deviations  from  an  example  more  nearly  repre¬ 
sentative  of  a  good  design,  the  velocity  at  the  blade  tip  for  example  II 
was  not  maintained  the  same  as  in  example  I,  and  consequently  the  differ¬ 
ences  between  examples  I  and  II  are  not  ascribable  solely  to  inlet  pre¬ 
rotation.  Comparison  of  the  effect  on  rotative  speed  (for  the  same 
turning  at  the  tip  and  the  same  work  output)  and  on  diffusion  at  the 
blade  roots  is,  however,  valid. 

The  results  of  the  computations  are  shown  in  figure  5.  In  this 
example,  the  turning  rate  is  predetermined  in  the  region  -0.5<z/rt<0 

by  the  velocity  distribution  upstream  of  the  leading  edge,  the  condition 
of  radial  blade  elements,  and  the  slope  of  the  leading  edge.  The 
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turning  rate  is  continuous  at  z  =  0,  and  satisfies  the  equation 

(0<z/rt  <l/3) 
(l/3<z/rt<l) 

Compared  with  the  first  example,  the  work- input  rate  at  the  blade  tip 
is  greater  in  the  initial  section;  however,  the  work  distribution  is 
still  unsatisfactory.  The  exit  blade  height  is  also  much  less,  with 
the  result  that  the  work  input  is  much  more  nearly  equalized;  at  the 
root,  the  value  of  AH/a^  is  1.8,  whereas,  at  the  tip,  the  value  is 
2.5.  These  values  correspond  to  stagnation  pressure  ratios  of  5.2  at 
the  root  and  8.35  at  the  tip.  This  variation  could  be  reduced  by 
sweeping  back  the  trailing  edge  to  provide  space  for  additional  work 
input  at  the  hub  streamline.  This  sweepback  is  more  feasible  with  the 
smaller  blade  height  and  work  variation  of  the  second  example  than  for 
the  first  example. 


_d  /rd£\ 
dz  \dz  / 


3.008  -  2.298  (3z  -  l) 


=  3.008 


[l  -  1/4  (3z  -  l)2] 


Example  III:  Effect  of  turning  rate.  -  In  order  to  reduce  the 
size  of  the  impeller  by  increasing  the  turning  rate,  example  III  was 
computed  with  twice  the  blade  turning  rate  of  example  II  for  z>0 
and  with  the  same  over- all  work  output.  One  result  is  the  greater 
diffusion  on  the  hub,  where  a  minimum  Mach  number  of  1.2  is  reached 
instead  of  1.4  as  in  the  previous  example  (See  fig.  6).  The  increased 
diffusion  probably  results  from  the  increase  in  hub  curvature,  which 
in  turn  is  caused  by  the  shorter  axial  distances  required  for  the 
changes.  Another  result  is  the  increased  height  of  the  blade  at  the 
exit,  and  an  increase  in  the  exit  Mach  number  from  2.3  to  2.7.  The 
increased  streamline  curvature  also  accounts  for  the  higher  exit 
velocity  and  this  higher  velocity  in  turn  explains  the  larger  blade 
height.  The  characteristic  line  through  the  casing  at  the  exit 
(r  =  r^,  z  =s  0.5  r^_)  intersects  the  hub  at  z  =  0.155  r^_.  The  char¬ 
acteristics  through  the  remaining  portion  of  the  hub  (z  >0.155  r^) 
are  affected  by  conditions  at  the  exit  (z  =  0.5  r^)  where  the  veloc¬ 
ities  increase  with  decreasing  radius.  The  result  is  a  large  increase 
in  velocity  at  the  blade  hub  in  the  region  z >0.155  r^. 

Example  IV:  Effect  of .diffusion  at  casing.  -  In  the  next  example, 
all  conditions  remained  the  same  as  in  the  third  example  except  that 
the  prescribed  velocity  at  the  .case  was  varied  smoothly  from  the 
value  of  1.94a  at  z  =  0.25  r^  to  1.40a  at  z  -  0.5  r^_.  The  desired 
reduction  has  been  obtained  in  blade  height  and  work-output  variation, 
but  the  discharge  velocity  has  a  large  value  of  2.6a  (fig.  7).  It  may 
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therefore  he  concluded  that  the  alternative  of  directly  prescribing  a 
condition  with  low  relative  velocities  at  the  exit  would  have  been 
more  effective  than  prescribing  conditions  at  the  blade  tip  for 
reducing  the  blade  height  at  the  exit. 

A  new  phenomenon  has  also  arisen  in  that  the  characteristics 
passing  through  z  =  0.5  rt,  r  =  0.92  rt;  and  z  =  0.5  r^,  r  =  0.893  rt 

intersect  at  z  =  0.397  rt,  r  =  0.832  rt>  At  this  point,  the  Mach 

number  changes  discontinuously  from  2.13  to  2.43  and  the  streamline 
inclination  to  the  axis  changes  from  14.9°  to  7.4°.  This  single  point 
does  not  represent  the  full  extent  of  this  difficulty.  Actually,  there 
is  a  region  bounded  by  a  curve  in  which  two  solutions  are  obtained; 
each  solution  is  determined  by  a  different  region  of  the  boundary.  Con¬ 
sequently,  the  boundary  conditions  are  incompatable  with  obtaining  a 
single-valued  continuous  solution.  This  situation  is  analogous  to  the 
Prandtl-Meyer  flow,  where  the  characteristics  intersect  at  the  corner. 

If  a  computation  were  started  for  a  simple  two-dimensional  supersonic 
flow  from  a  streamline  shape  and  from  the  velocity  distributed  thereon, 
which  were  exactly  the  same  as  on  one  of  the  streamlines  of  a  Prandtl- 
Meyer  flow  around  a  corner,  then,  when  the  solution  was  extended  toward 
the  corner,  the  characteristics  would  intersect  at  the  corner  in  the 
manner  of  figure  7.  The  intersection  on  figure  7  can  be  taken  to  repre¬ 
sent  a  physical  boundary  to  the  fluid  in  that  the  streamline  passing 
through  that  point  would  be  taken  for  the  hub  shape.  Only  90  percent 
of  the  flow  (ty  =  0.9)  could  be  accommodated  in  the  impeller  of 
example  IV  if  continuous  flow  with  the  imposed  boundary  conditions  were 
required.  An  obvious  method  of  increasing  the  flow  is  to  increase  the 
slope  of  the  characteristic  lines  at  the  exit  for  r<0.9  r.  .  This 

increase  could  be  accomplished  by  decreasing  the  relative  velocities 
there  and  discarding  the  condition  v  =  0  to  obtain  a  shorter  blade 
height . 

Example  V:  Increasing  blade  thickness.  -  In  none  of  the  examples 
computed  so  far  was  blade  thickness  considered.  An  obvious  effect  of 
increasing  blade  thickness  is  to  decrease  the  velocity  in  a  region 
where  there  is  no  independent  control  of  the  velocity  (such  as  the 
region  z<0.35  rt  in  example  I,  fig.  4).  This  phenomenon  provides 
the  degree  of  freedom  required  in  such  regions  to  control  the  relative 
velocity  |W|  as  well  as  the  surface  velocities  |Wg|  and  |Wp|  .  In 

regions  where  the  velocity  |W|  at  the  blade  tips  is  under  control 
(the  entire  tip  sections  for  examples  II,  III,  IV) ,  altogether  dif¬ 
ferent  effects  may  be  expected.  The  effect  of  thickness  under  these 
circumstances  was  investigated  by  computing  another  example  with  a 
cosine-like  bump  in  the  blade  thickness  starting  at  z  =  0.25  rt . 
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(See  the  previously  presented  table  of  boundary  conditions.)  The  maxi¬ 
mum  thickness  was  reached  at  z  =  0.375  r^ .  At  r  =  r^  the  thickness 

is  20  percent  of  the  blade  spacing,  and  at  the  hub  the  value  is  25  per¬ 
cent.  The  results  of  the  computation  are  shown  in  figure  8.  The  effect 
of  the  bump  is  not  felt  upstream  of  the  characteristic  with  positive 
slope  passing  through  z  =  0.25  rt,  r  =  r^.  The  most  obvious  effect  is 

a  decrease  of  flow  Mach  number  in  the  region  of  increasing  thickness 
(and  inward  from  the  casing,  where  the  velocity  is  not  preassigned) ,  as 
might  be  expected  from  one-dimensional  theory.  At  z  =  0.25  r^_  and 

z  =  0.375  r^,  the  change  in  blade  height  is  about  10  percent  as  compared 

with  example  IV.  Another  effect  of  thickness  is  to  be  seen  in  the 
waviness  of  the  hub  streamline  in  the  region  0.15<  z/r^c  0.25  where 

the  thickness  has  not  yet  started  to  increase.  This  bump  may  be  traced 
in  other  streamlines  along  the  characteristic  running  through  the  point 
of  maximum  thickness  at  the  case  (r  =  r^,  z  =  0.375  r^_) .  The  effect 

of  the  bump  on  the  hub  is  therefore  to  increase  the  local  velocity  by 
curvature  in  order  to  hold  the  required  velocity  at  the  casing;  this 
action  from  hub  to  casing  takes  place  along  the  characteristic  line  as 
simple  waves  in  two  dimensions,  and  the  velocity  is  seen  to  have 
increased  all  along  this  characteristic  line.  The  intersection  of  the 
characteristics  and  the  double-valued  solution  occur  near  the  exit  over 
a  considerably  larger  portion  of  the  flow  than  in  the  previous  example. 
In  this  example  also,  if  a  lower  velocity  were  prescribed  at  the  exit, 
a  unique  solution  would  be  obtained.  It  would  also  be  feasible  to  carry 
out  this  idea  and  at  the  same  time  to  prescribe  another  condition  (such 
as  v  =  0)  provided  that  one  were  willing  to  lose  control  of  the  condi¬ 
tion  at  the  back  end  of  the  casing  (r  =  r^)  in  a  manner  analogous  to 
that  used  at  the  front  end  of  example  I. 


Example  VI:  Effect  of  leading-edge  sweep.  -  It  is  of  more  technical 
interest,  however,  to  carry  out  this  scheme  at  the  entrance,  for  then 
the  leading-edge  sweep  may  be  reduced  from  50  percent  of  the  axial  depth 
of  the  machine  to  considerably  less.  A  computation  of  the  inlet  flow 
conditions  for  such  an  impeller  is  shown  in  figure  9,  where  the  condi¬ 
tions  at  the  leading  edge  are  the  same  as  for  example  III.  The  blade 


shape  is  different,  with  the  second  derivative  of  angle 


a2e  _  a2f(z) 

.  2  “  “  .2 

dz  dz 

constant.  Also,  the  velocity  at  r  =  1.0  cannot  be  prescribed.  Only 
in  the  region  enclosed  by  the  casing  and  leading- edge  characteristics 
was  the  computation  completed.  The  result  (fig.  9)  shows  a  large 
increase  in  Mach  number,  especially  at  the  blade  tip.  Because  of  this 
increase,  the  blade  turning  accomplishes  no  work  in  the  beginning  of 
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the  impeller,  but  shows  a  region  of  energy  abstraction  from  the  air;  the 
Initial  negative  rotational  component  becomes  more  negative  as  the  air 
turns  axially  and  the  velocity  increases  as  in  example  I.  The  situation 
is  quite  different  at  the  back  end  where  the  blade  angles  approach  zero 
and  the  tangential  component  does  likewise;  all  the  work  input  is  con¬ 
centrated  in  the  back  end  of  this  impeller  as  in  the  other  examples. 
Postponement  of  work  input  to  the  back  end  of  the  impeller  is  to  be 
expected  in  any  case  where  the  velocity  is  increasing  as  the  blades 
turn  toward  the  axial  direction,  or  if  the  velocity  decreases  as  the 
blades  turn  past  the  axial  direction.  Internal  flow  and  work  distribu¬ 
tion  could  be  considerably  improved  in  this  example  by  increasing  blade 
thickness  in  the  region  of  large  velocity  increase.  * 


CONCLUSIONS 

General  equations  developed  for  isentropic,  axisymmetric  flows  in 
rotating  impellers  with  any  given  blade  shapes  showed  that  the  blade- 
force  term  can  be  eliminated  and  solved  for,  if  desired,  after  the  flow 
solution  is  completed.  Further,  if  the  flow  upstream  of  the  impeller 
is  vortex- free,  then  a  potential  function  of  only  radial  and  axial 
distances  exists,  which  may  be  used  to  find  the  projection  of  the  abso¬ 
lute  velocity  on  the  blade  surface. 

For  supersonic  flow,  the  method  of  characteristics  was  used  to 
compute  the  flow  for  several  examples  with  radial  blade  elements  and 
vortex-free  inflow.  These  examples  show  that 

(1)  With  zero  radial  inflow,  the  influence  of  the  leading  edge  and 
the  blade  turning  on  the  casing  velocity  in  regions  where  this  velocity 
cannot  be  prescribed  in  advance  causes  an  increase  in  velocity  and  a 
reduction  of  blade  loading  and  work  input  in  the  initial  portion  of  the 
compressor.  These  effects  can  be  compensated  for  to  some  extent  by 
increasing  blade  thickness. 

(2)  Sufficient  leading-edge  sweepback  adds  a  degree  of  freedom  to 
the  boundary  conditions  at  the  blade  tips,  and  therefore  permits  a  more 
constant  or  a  decreasing  velocity  inside  the  impeller. 

(3)  Reduction  in  the  exit  velocity  results  in  a  shorter  blade 
height  and  a  more  uniform  distribution  of  work  input  at  the  exit. 

(4)  Increase  in  the  rate  of  turning  causes  more  rapid  hub  curva¬ 
tures  and  consequently  larger  internal  velocity  gradients.  Diffusion 
on  the  hub  is  therefore  greater  in  the  initial  part  of  the  impeller, 
and  a  lower  velocity  minimum  is  reached  there,  whereas  the  gas  velocity 
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at  the  huh  near  the  exit  reaches  larger  values.  This  situation  requires 
larger  hlade  heights  at  the  exit  with  a  larger  exit  velocity  at  the  root 
and  greater  variation  of  work  output  from  root  to  tip.  Diffusion  of  the 
gas  velocity  at  the  tip  will  moderate  this  effect,  as  will  also  a  relax¬ 
ation  of  the  condition  zero  radial  velocity  at  the  exit  and  the  substi¬ 
tution  of  a  condition  giving  lower  exit  velocities. 

(5)  Not  every  set  of  boundary  conditions  results  in  a  continuous, 
single-valued  solution  in  the  region  of  interest.  In  the  examples  con¬ 
sidered,  a  relaxation  of  the  condition  zero  radial  velocity  at  the  exit 
would  have  permitted  another  condition  giving  lower  velocities  there. 
This  condition  wodld  have  eliminated  the  region  of  the  double-valued 
solution  and  given  a  shorter  blade  and  a  more  nearly  equalized  work 
output . 

(6)  The  effect  of  variations  in  blade  thickness  in  regions  where 
velocity  is  also  prescribed  at  the  blade  tip  is  to  change  the  velocity 
as  in  one- dimensional  flow  but  not  by  the  same  amount,  for  there  is 
some  compensation  by  variations  in  blade  height.  The  prescribed  veloc¬ 
ity  and  thickness  require  hub  curvatures  that  cause  velocity  changes  to 
be  projected  like  simple  waves  on  characteristic  curves.  It  is  there¬ 
fore  possible  to  find  regions  of  increased  velocity  when  the  blade 
thickness  is  increased. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  January  24,  1951. 
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APPENDIX  -  SYMBOLS 


The  following  symbols  are  used  in  this  report: 
r  -  Nt/2?t 

sonic  velocity  of  gas 

sonic  velocity  of  gas  at  conjunction  of  leading  edge  and 
impeller  case 

critical  speed  of  gas  upstream  of  impeller 
an  abbreviation  (see  equations  (21a)) 
abbreviations  (see  equations  (21a)) 


sum  of  stagnation  enthalpy  and  potential  energy  of  gas  (specific) 

/  w2  2  21 

in  moving  coordinates  [=h  +  —  - 

\  Li 


-) 


distributed  blade  force  per  unit  mass  of  fluid  (a  vector) 

part  of  blade  surface  function  dependent  on  r  and  z  onlv 
(a  =  6  +  f(r,z)) 

an  abbreviation  (see  equations  (21a)) 

p 

absolute  stagnation  enthalpy  (h  +  V  /2) 

work  input  to  gas  (rise  in  absolute  stagnation  enthalpy) 

static  enthalpy  of  gas  (specific) 

unit  vector  parallel  to  curve  9  =  constant,  z  =  constant 
(i  =  7r) 

an  abbreviation  (see  equations  (21a)) 

unit  vector  parallel  to  curve  z  =  constant,  r  =  constant 

(j  =  rye) 

unit  vector  parallel  to  curve 
(k  =  yz) 


k 


r  ss  constant,  9  =  constant 
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L  an  abbreviation  (see  equations  (21a)) 

M  local  relative  Mach  number 

N  number  of  blades 

p  static  pressure  of  gas 

q  unit  vector  parallel  to  W 
R  vector  of  r 

r  normal  distance  from  axis  of  rotation 

r,  distance  from  axis  of  rotation  to  blade  tip  at  leading  edge 

U 

ds  arc- length  element  parallel  to  W 

t  blade  thickness  measured  in  direction  of  rotation 

u  rotational  component  of  relative  gas  velocity 

V  absolute  gas  velocity  (a  vector)  |iv  +  j(u-Knr)  +  kw] 

v  radial  component  of  gas  velocity 

W  relative  gas  velocity  (a  vector)  (iv  +  ju  +  kw) 

w  axial  component  of  gas  velocity 

z  distance  measured  parallel  to  axis  of  rotation 

a  function  for  blade  surface  (a(r,0,z)  =  constant) 

p  angle  between  meridional  projection  of  leading  edge  and  axis  of 
rotation 

X  ratio  of  specific  heats 

f;  either  or  in  general  form  of  equations 

£+  slope  of  curve  for  |  variable,  rj  constant 
slope  of  curve  for  rj  variable,  i  constant 
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T)  one  of  characteristic  coordinates  for  gas  flow  in  impeller 
9  meridional  angle  (for  cylindrical  coordinates) 

H  absolute  moment  of  momentum  (r(u  +  or)) 

t  one  of  characteristic  coordinates  for  gas  flow  in  impeller 
p  gas  density 

a  either  tj  or  %  in  general  form  of  equations 

$  velocity-potential  function  of  r  and  z  only  that  exists  for 
initially  irrotational  absolute  motion 

\|r  stream  function  (yt|/  =  p  AjxV) 

Q  vector  form  of  cjd 

cn  angular  velocity  of  impeller 

Subscripts : 

0  upstream  of  impeller 

p  pressure  surface  of  blade  at  tip  . 

s  suction  surface  of  blade  at  tip 

Subscripts  r,  0,  z,  o,  T|  indicate  partial  differentiation. 
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Figure  2.  -'Relation  between  shape  of  boundary  and  boundary  conditions 
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Center  line 

(a)  Conditions  on  b’B  and  BC  unrelated. 


(b)  Conditions  on  B’B  related  to  those  on  BC . 


Figure  3.  -  Boundary  conditions  for  impeller. 
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